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■ The mass of the top quark has been measured at Fermilab using the Tevatron pp collider to a 

£N| precision of less than 1%. I discuss the individual measurements that go into this result and 

its impact on our understanding of the standard model of the electroweak interactions. 
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The top quark is the most massive fundamental particle that we know today. It is the 
dominant contributor to radiative corrections for many standard model processes. In order to 
make precise predictions based on the standard model we have to therefore know its mass to 
high precision. A well-known example for this is the W boson mass. The dominant radiative 
corrections to the W boson mass come from loops containing top and bottom quarks, which are 
proportional to the top quark mass squared, and from loops containing Higgs bosons, which are 
logarithmically dependent on the Higgs boson mass. Thus for every value of the Higgs boson 
mass rrtff, the standard model defines a unique line in the plane defined by W boson mass and 
top quark mass shown in Figure [2 Vice versa, if we measure the W boson and top quark masses 
precisely we can constrain the allowed parameter space for the standard model in this plane and 
therefore limit the allowed range for the Higgs boson mass. 

We study top quarks at the Fermilab Tevatron pp collider at \fs = 1.96 TeV. The dominant 
production mechanism for top quarks is top-antitop quark pair productions via quark-antiquark 
annihilation or gluon fusion. Within the standard model the top quark decays almost exclusively 
to a W boson and a b quark. The decay modes of the W boson then define the final state 
signatures. In 6% of all ti decays both W bosons produce electrons or muons in their decay - 
we call this the dilepton channel. In 38% of all ti decays one W boson produces an electron or 
a muon in its decay - this is the lepton+jets channel. The remaining 56% of ti decays contain 
no electrons or muons from W decays - we call this the all-jets channel. In addition to the W 
bosons, all ti decays also produce b and b quarks that give rise to jets. Since 6-flavored hadrons 
are long-lived, we can identify such jets by the presence of tracks from their secondary decay 
vertices that do not point back to the primary interaction point. Jets with such tracks or with 
explicitly reconstructed secondary vertices likely originate from the fragmentation of a b quark 



and are called 6-tagged. Requiring at least one 6-tagged jet in the event is a good way to increase 
the signal-to-background in the data. 

The measurement of the top quark mass relies on the precise reconstruction of the energy of 
its decay products. We calibrate the measurement of the energies of electron and muons using 
Z — ► ee or Z — > ///i decays. The calibration of the quark energies is a bigger challenge. Quarks 
manifest themselves as jets in the detector and the first step is to calibrate the measurement of 
the energy in the jet cone to the energy of all particles from the hadronization of the original 
quark in the same cone. The response of the calorimeter to the particles is calibrated by 
balancing the transverse momentum in events with a direct photon and a jet. Then Monte 
Carlo-derived corrections are applied to take this calibration to the parton level. Dijet events 
are used to make the calorimeter response to jets uniform in the calorimeter at all momenta. We 
call this the external jet energy calibration and it is sytematically limited at a precision of about 
1.7%. Alternatively, we can use hadronic decays of W bosons from top decay to calibrate the 
jet energy scale with top decay events. This determines an additional overall calibration factor 
for jet energies. We call this the in situ calibration. This calibration is presently statistically 
limited by the number of hadronic W decays. 

Two basic types of algorithms are currently used for the top mass measurements. The 
template methods reconstruct one or more observables from each event that serve as estimators 
for the top quark mass and then compare their distributions to Monte Carlo-derived distributions 
for a range of top quark masses, the so called templates. More recently, methods have been 
developed that derive a likelihood as a function of the top quark mass for each event. This has 
the advantage that every event can be weighted according to its resolution and probability to be 
from tt decays. The most prominant example for the latter is the matrix element method that 
computes the probability for an event to be produced in tt decays as a function of the top quark 
mass based on inegrating over the tree-level matrix element for tt production and decay^l. 

It has become standard procedure to calibrate the performance of the measurement al- 
gorithms using Monte Carlo generated pseudoexperiments. Events are drawn randomly from 
simulated signal and background samples according to the number of events observed in the 
data and the expected or observed signal-to-background. Each pseudoexperiment is analyzed 
with the same algorithm as the data. The average and root mean square of the measured masses 
and their pull distribution characterize the performance of the algorithm. 

Dilepton decays contain two 6-jets, two charged leptons and several neutrinos in the final 
state. This channel has the advantage of a high signal-to-background (3:1) even without using b- 
tagging. Its disadvantage is that the number of measured quantities in these events is insufficient 
to constrain the kinematics of the events. 

DO uses two different algorithms to measure the top quark mass in dilepton events based on 

1 fb _1 of data. The matrix weighting algorithm finds for every top mass the possible top and 
antitop momenta and then assigns a weight to the events based on parton distribution functions 
and the lepton energy distribution. The value of the top quark mass at which this weight is 
maximized is used as the estimator. Its distribution is fit to templates to extract a top quark 
mass of 175.2±6.1(stat)±3.4(syst) GeV. The neutrino weighting algorithm loops over a range 
of pseudorapidities for the two neutrinos and uses a \ 2 variable based on the consistency of 
the sum of the neutrino momenta with the observed missing p? as weight. A two dimensional 
template fit to the distributions of mean and rms of the weight curves gives a top quark mass 
of 172.5±5.8(stat)±3.5(syst) GeV. Combining the two gives 173.7±5.4(stat)±3.4(syst) GeV^. 

CDF use the matrix element method in the dilepton channel. The analysis is carried out on 

2 fb _1 of data and optimizes the selection with an evolutionary neural network to minimize the 
statistical uncertainty. The result is 171.2±2.7(stat)±2.9(syst) GeV^. CDF also apply the neu- 
trino weighting method. Based on 1.9 fb" 1 and a two-dimensional template fit to the mass with 
the largest weight and to the scalar sum of the jet pts gives 171.6±3.3(stat)±3.8(syst) GeV^. 



In lepton+jets decays we measure everything except the component of the neutrino momen- 
tum along the beam direction. With the two mass constraints on the decay products of the W 
bosons and the requirement for both top quark masses to be equal, a kinematic fit with two 
constraints is possible. Signal-to-background is typically 1:2 without 6-tagging and 2:1 after 
requiring at least one 6-tagged jet. 

DO use the matrix element method to measure the mass in the lepton+jets channel. DO 
integrate over the tree-level matrix elemnts for tt-production and for VF+jets production to 
compute the likelihood that the event originates from either signal or background. The jets 
from the hadronic W boson decay are used to constrain the jet energy scale calibration. The 
result rat = 170.5 ± 2.5(statffijes) ± 1.4(syst) GeV from the first 0.9 fb _1 of data from Run II 
has already been presented at conferences last year. DO have analyzed an additional 1.2 fb _1 of 
data resulting in nit = 173.0 + 1.9(stat©jes) ±1.0(syst) GeV. DO combine the two measurements 
to obtain m t = 172.2 ± l.l(stat) ± 1.6(syst©jes) GeV ? . 

CDF also use the matrix element method with in situ calibration of the jet energy scale. 
CDF only integrate over the ti-matrix element to compute the likelihood to observe the event 
from it-production. No integration over the background matrix elements is performed. CDF 
use a neural network discriminant to determine an event-by-event signal-to-background ratio 
and then subtract the bias in the likelihood expected from the background events in the sample. 
Based on 1.9 fb _1 CDF measure m t = 172.7±1.8(stat©jes) ±1.2(syst) GeV^. Alternatively CDF 
also perform a template fit to the best top quark mass and the best dijet mass from a kinematic 
fit of the event to the tt- hypothesis on the same data set. The dijet mass constrains thejet 
energy scale. With this analysis CDF measure mt = 171.8 ± 1.9(stat©jes) ± l.O(syst) GeV^H 

The all jets channel has the largest branching fraction but also by far the largest background. 
CDF use events with 6-8 jets, of which at least one has to be 6-tagged. They fit the best top 
quark mass and dijet masses from a kinematic fit of the leading six jets to the tt- hypothesis to 
templates and obtain m t = 177.0 ± 3.7(stat©jes) ± 1.6(syst) GevE. 

We use the BLUE method^ to combine the best top quark mass measurements from Runs I 
and II in each channel by the to experiments, including the correlations between the systematics 
uncertainties of these measurements. Fig. [T] compares these measurements and Table [U summa- 
rizes the Run II measurements that enter in this combination and their systematic uncertainties. 
It can be seen that statistical and systematic uncertainties are of similar magnitude and that the 
largest source of systematic uncertainty is the jet energy scale calibration. The world average 
top quark mass is mt = 172.6 ± 1.4 GeV, which is for the first time below 1% in precision. The 
top quark masses determined in the three different final state channels are 169.8 ±3.1 GeV in 
the dilepton channel, 172.4 ± 1.5 GeV in the lepton+ jets channel, and 177.3 ± 3.9 GeV in the 
all-jets channel, all consistent with the overall average^l 

The global electroweak fit was repeated with this new world average for the top quark mass 
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as input 1 ^. The x is minimized at a value of 17.2 for 13 degrees of freedom, which corresponds 

to a probability of 19% for a Higgs boson mass of 87±f 7 GeV. The 95% C.L. upper limit for 

the Higgs boson mass is 160 GeV. Figure [2] shows the measured values of the masses of the top 

quark and the W boson^l. 

In conclusion, we have measured the top quark mass with 0.8% precision to be 172.6+1.4 GeV 
using about 2 fb _1 of data from the Fermilab Tevatron collider. By the end of Run II we expect 
to accumulate a total of around 6-8 fb _1 which will allow us to push the experimental precision 
of this measurement below 1 GeV. 



I would like to thank the organizers for inviting me to present these results at this very 
exciting and stimulating conference. 
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Figure 1: Measurements of the top quark mass. 



Figure 2: W boson mass versus top quark mass 
plane. The shaded band is consistent with the stan- 
dard model for 114 GeV< m H < 1 TeV. 



Table 1: Run II top quark mass measurements and systematic uncertainties used in the world average. 



channel 


dilepton 


1+jets 


all-jets 


T LUJ 


world 




experiment 


CDF 


DO 


CDF 


DO 


CDF 


CDF 


average 




JCdt 


2.0 


1.1 


1.9 


2.1 


1.9 


0.7 




fb- 1 


result 


171.2 


173.7 


172.7 


172.2 


177.0 


180.7 


172.6 


GeV 


jet energy calibration 


1.5 


3.1 


1.5 


1.3 


2.0 


0.3 


0.9 


GeV 


signal model 


0.7 


0.8 


0.6 


0.7 


0.6 


1.4 


0.5 


GeV 


background model 


0.4 


0.6 


0.6 


0.4 


1.0 


7.2 


0.4 


GeV 


fit 


0.6 


0.9 


0.2 


0.1 


0.6 


4.2 


0.1 


GeV 


Monte Carlo model 


0.7 


0.2 


0.4 


0.0 


0.3 


0.7 


0.2 


GeV 


syst uncertainty 


2.8 


3.4 


1.7 


1.6 


2.4 


8.5 


1.1 


GeV 


stat uncertainty 


2.7 


5.4 


1.2 


1.1 


3.3 


14.5 


0.8 


GeV 


total uncertainty 


3.9 


6.4 


2.1 


1.9 


4.1 


16.8 


1.4 


GeV 
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